FLORENTIN SMARANDACHE 
Another Integer Algorithm 
To Solve Linear Equations 
(Using Congruences) 


In Florentin Smarandache: “Collected Papers”, vol. I (second edition). Ann 
Arbor (USA): InfoLearnQuest, 2007. 


ANOTHER INTEGER ALGORITHM TO SOLVE LINEAR EQUATIONS 
(USING CONGRUENCES) 


In this section is presented a new integer number algorithm for linear equation. 
This algorithm is more “rapid” than W. Sierpinski's presented in [1] in the sense that it 
reaches the general solution after a smaller number of iterations. Its correctness will be 
thoroughly demonstrated. 


Another Integer Algorithm. 
Let's us consider the equation (1); (the case a,b €Q, i=1,n is reduced to the 


case (1) by reducing to the same denominator and eliminating the denominators). Let 
d=(a,,...,a,). If d Ib then the equation does not have integer solutions, while if d |b 
the equation has integer solutions (according to a well-known theorem from the number 
theory). 
If the equation has solutions and dz we divide the equation byd. Then d=1 
(we do not make any restriction if we consider the maximal co-divisor positive). 
Also, 
(a) If all a, the equation is trivial; it has the general integer solution 
x, =k, €Z, i=1,n, when b=0 (the only case when the general solution is 
n-times undetermined) and does not have solution whenb # 0. 
(b) If ai, 1<i<n such that a,=+1 then the general integer solution is: 


ie \ 
x, na] Za -b and x, =k, €Z, se{I,...,n}\ {i} 


The proof of this assertion was given in [4]. All these cases are trivial, therefore 
we will leave them aside. The following algorithm can be written: 


Input 
A linear equation: 


Oy “Xa babe apela 
i=l 


with not all a, = 0 and (a,,...„a,)=1. 


Output 

The integer general solution of the equation. 
Method 

l. hkh=1,p=1 


r 











r| < a) and determine r and the pair 


2. Calculate min | ; 
1<i,j<n 


, r=a,(moda,), 


(i, j) for which this minimum can be obtained (when there are more possibilities we have 
to choose one of them). 
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3. If |r| # go to step 4. 
If |r| =1, then 








(A) Substitute the values thus determined of these unknowns in all the 
statements (p), p = 1,2,... (if possible). 


(B) From the last relation (p) obtained in the algorithm substitute in all 
relations: (p —1),(p —2),...,(1) 

(C) Every statement, starting in order from (p = 1) should be applied the 
same procedure as in (B): then (p — 2), ...„(3) respectively. 

(D) Write the values of the unknowns x, i= In, from the initial 


equation (writing the corresponding integer number parameters from 
the right term of these unknowns with k,,...,k,_, ), STOP. 


i 





4. Write the statement (p): x, =t,- x 
j 
5. Assign x, =t, h=h+1 
a, =r p=ptl 


The other coefficients and variables remain unchanged go back to step 2. 
The Correctness of the Algorithm 


Let us consider linear equation (2). Under these conditions, the following 
properties exist: 


Lemma 1. The set M =| \r 


Proof: 

Obviously M c N' and M is finite because the equation has a finite number of 
coefficients: n, and considering all the possible combinations of these, by twos, there is 
the maximum AR? (arranged with repetition) = n” elements. 

Let us show, by reductio ad absurdum, that M +Ø. 

Mz0 & a,=O(moda,) Vi, j = ln. Hence a, = 0(moda,) Vi, j = Lin. Or this 





„r=a,(moda,), 0<|r| < a, has a minimum. 


is possible only when a, „Vi, j=Ln, which is equivalent to 
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(a,,...4,)=a,, Vi eLn. But (a,,..,a,)=1 are a restriction from the assumption. It 
follows that la, = Ln, Vi el,n a fact which contradicts the other restrictions of the 


assumption. 
M +0 and finite, it follows that M has a minimum. 


„Vieln. 








Lemma 2. If |r| = min M , then Ir] <la, 
Si,j<n 
Proof: 


We assume conversely, that 3i,, 1 <i, <n such that |r| 2 








a. 
10 


Then | > min | a, = a, 


1<j<n 





#1, 1<j,<n.Leta,,1< p, <n, such that la, |> la, and 
. ee 0 
a,, is not divided by a,. 


There is a coefficient in the equation, 





d; | which is the minimum and the coefficients 
are not equal among themselves (conversely, it would mean that (a,,..,a,)=a,=+1 
which is against the hypothesis and, again, of the coefficients whose absolute value is 
greater that la,, | not all can be divided by a, (conversely, it would similarly result in 


(a,,...4,) = a, # +1. 
We write |a, /a,, ] =q, EZ (integer portion), and r=a, —qa,, €Z. We have 


a,|<|r|. Thus, we have found an 7, which 


Po 10 








a, = r(moda,) and 0< Inl < a,,| < 
In| <|r| contradicts the definition of minimum given to |r]. 


Thus |r| <|a, „Vi el,n. 





Lemma 3. If |r| = min M = 1 for the pair of indices (i, j), then: 








x, =k, eZ, se {1,...,2}\{i, j} 


is the general integer solution of equation (2). 
Proof: 


95 


Let 22a, e=],n, be a particular integer solution of equation (2). Then 














Jk, =x) €Z, s Ee {,....a}\ {i,j} and t, = x94 Si e EZ (because a, —r = Ma, ) such 
a; 
that: 
| 0 iTr | = 0 0 
x, =r-a;| x; + A = a,x, +tb=x 
a; s=l 
sefi, j} 
0 ia 0 a —Lr n a; 0 
norma [ats a} > a,x, + b=x, 
a, a, T a, 


and x, =k, =x), s€ {1,.. n} \ {i j}. 


Lemma 4. Let |r| # and (i, j) be the pair of indices for which this minimum can 
be obtained. Again, let’s consider the system of linear equations: 


n 
at, rX; + y ax, =b 


(3) sefi, j} 





Then x, =x?, e= l,n is a particular integer solution for (2) if and only if Ce ae 


a; — 





Pas . : . 
x; is the particular integer solution of (3). 


l 


ee{l,.. n} \{j} and t, =t =x; + 


v 


Proof: 
x, =x?, e=1,n is a particular solution for (2) if and only if 
n n a. SI r 
Sax =0S $ ax +a, x + — x +r =b 
e=l 


a; 





s=l 
s¢{i,j} 


0 r 


n 
es 
O ad, HIX, + > a,x? =b and t, =x] + 
s=1 j 


j 
sefi, j} 





0 
x, EZ Sx, =%,5 


ee {I,...,.n}\{j} and 1, =t; is a particular integer solution for (3). 


Lemma 5. The previous algorithm is finite. 

Proof: 

When |r| =1 the algorithm stops at step 3. We will discuss the case when 
r| eN'”. We will show that the row of 


r—s successively obtained by following the algorithm several times is 
decreasing with cycle, and each cycle is not equal to the previous, by 1. Let 7, be 


|r| #1. According to the definition of r, 
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the first obtained by following the algorithm one time. Iz #1 then go to step 4, 
, Vi=l,n. 


Now we shall follow the algorithm a second time, but this time for an 
equation in which 1 (according to step 5) is substituted by a,. Again, according 


and then step 5. According to lemma 2, a; 





n| < 








to lemma 2, the new |r| written |r,| will have the propriety: |r,| < |7]. We will get 
to |r| = 1 because |r|>1 and |r|<oo, and if |r| #1, following the algorithm once 
again we get |r|< Izl and so on. Hence, the algorithm has a finite number of 
repetitions. 


Theorem of Correctness. The previous algorithm calculates the general 
solution of the linear equation correctly (2). 

Proof: 

According to lemma 5 the algorithm is finite. From lemma 1 it follows 
that the set M has a minimum, hence step 2 of the algorithm has meaning. When 
|r| = 1 it was shown in lemma 3 that step 3 of the algorithm calculates the general 


integer solution of the respective equation correctly the equation that appears at 
step 3). In lemma 4 it is shown that if |r|#1the substitutions steps 4 and 5 


introduced in the initial equation, the general integer solution remains unchanged. 
That is, we pass from the initial equation to a linear system having the same 
general solution as the initial equation. The variable h is a counter of the newly 
introduced variables, which are used to successively decompose the system in 
systems of two linear equations. The variable p is a counter of the substitutions of 
variables (the relations, at a given moment between certain variables). 

When the initial equation was decomposed to |r| = 1, we had to proceed in 


the reverse way, i.e. to compose its general integer solution. This reverse way is 
directed by the sub-steps 3(A), 3(B) and 3(C). The sub-step 3(D) has only an 


aesthetic role, i.e., to have the general solution under the form: x, = f/(k,,...k,_1); 


i=l,n, f, being linear functions with integer number of coefficients. This “if 


possible” shows that substitutions are not always possible. But when they are we 
must make all possible substitutions. 


Note 1. The previous algorithm can be easily introduced into a computer 
program. 


Note 2. The previous algorithm is more “rapid” than that of W. 
Sierpinski’s [1], i.e., the general integer solution is reached after a smaller number 
of iterations (or, at least, the same) for any linear equation (2). 

In the first place, both methods aim at obtaining the coefficient +1 for at least 
one unknown variable. While Sierpinski started only by chance, decomposing the 
greatest coefficient in the module (writing it under the form of a sum between a 
multiple of the following smaller coefficient (in the module) and the rest), in our 
algorithm this decomposition is not accidental but always seeks the smallest |r| 
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and also choose the coefficients a, and a, for which this minimum is achieved. 


That 1s, we test from the beginning the shortest way to the general integer 
solution. Sierpinski does not attempt to find the shortest way; he knows that his 
method will take him to the general integer solution of the equation and is not 
interested in how long it will take. However, when an algorithm is introduced into 
a computer program it is preferable that the process time should be as short as 
possible. 


Example 1. 

Let us solve in Z? the equation 17x- 7y+10z=~-12. 
We apply the former algorithm. 

|. h=1,p=l 

2.r=3,i=3,j=2 

3. |3|#1 go on to step 4. 





10-3 
4. (1) y=4,- -Z=H +z 
5. Assign 
y =t h:=2 
a,=3 p=2 


with the other coefficients and variables remaining unchanged, go back to 
step 2. 
2. r=-l, i=l, j=3 
3. |-1]=1 
x=-l(-3t, - (7t )- 12) = 3t, — 7t — 12 
17-(-1) -1-17 
3 j 3 


(A) We substitute the values of x and z thus determined into the only 
statement (p) we have: 


(A) y=t,+2=--17t,4+43t,-72 








2-1 171, +( Tt): ( 12) =17t, + 42t, — 72 


(B) The substitution is not possible. 
(C) The substitution is not possible. 
(D) The general integer solution of the equation is: 


x = 3k, -7k, +12 
y = -I7k, + 43k, -72 
z=-17k, +42k,-72; kk, eZ 
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